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Abstract 



We describe an ansatz for symmetry reduction of the Lane-Emden 
!>■ I equation for an arbitrary polytropic index n, admitting only one sym- 

metry generator. For the reduced first order differential equation it 
{/^ . is found that standard reduction procedure do not admit any non- 

^D I trivial Lie point symmetry. However some special solutions for the 

t^^ ' differential equation are obtained. 



The simplest model of a gravitationally bound spherically symmetric 
polytrope with a power law distribution function is determined by the Lane- 
Emden equation [H [2] 

r"^ ar ar 

where ?/'(r) is the dimensionless gravitational potential, r is the dimensionless 
radius of the sphere, and n is the polytropic index. For the second order 
differential equation the symmetry generator has the form[3] 

Here the prime denotes differentiation with respect to the independent vari- 
able r. For the differential equation the symmetry condition is given by 

1 



^ojir^Tp^Tp') = rj"{r-^ip), raodip" = uj{r,ilj,i/j'). (3) 

This gives an identity in the dependent and independent variables {r^ip). 
From (1) we have 

uj{r,i),ij') = -ij{rY--ij'. (4) 

r 

The infinitesimals (^, rj) are now determined by the identity 



2 2 2 

We can determine the infinitesimals C, and r^ by an iterative procedure. First 
equating coefficients of {ip')^ and {ip'^ in (5), we find that ^ 7^ 'C('^) ^-^^d 
^,V'V ~ 0- Therefore we assume that C, = <y{r) and rj = j3{r)il) + 7(r). The 
symmetry condition then reduces to 



(^" + -^')[p - 2a'] + \ilj'a - mP^'-'i/Stlj + 7) - -[/^V + 7' + ^'(/5 - a')] 



= (3"^ + ^" + ^'[2(3' -a"] (6) 

Equating the coefficients of ip' and ■?/'" we obtain the following two conditions 

2a - 2ra' - [2/3' - a'Y = 0, [(3- 2a'] - nf3 = (7) 

The coefficients of r^ equated to zero leave for the first two terms in the first 
condition a = kr where k does not depend on r.The second equation gives 
(3 = 2a' /{I — n). Substitution back into the determining equation yields 

n7^"-^ + -7' + 7" = (8) 

r 

Comparing coefficients of ^'"'~ we find that 7 = 0. Thus the infinitesimals 
are 

—2k 

{^,r]) = {kr, -^), n^l. (9) 

n — 1 



We have here a one-parameter group Gi of scahng transformation with a 
parameter k. Thus reduction in the order of the equation by one is possible 
[1]. For convenience we take k to be unity. 

Let [t, s) be the canonical variables for the differential equation (1) with 
infinitesimals (9). Then 

s = Jj, t = j{idy-r,dx)- (10) 



gives 



2 
s = Inx, t = |/exp( -); n ^ 1. (11) 



We have for s = s{t) and s' = ds/dt: 



If 

dr = d{e') = e'ds, d^ = d{te^''^^-''^) = e'^'/'-^-'^Ut + e^'/^^'^'Us. 



1 — n 



(12) 



which implies that 



and 



dih e^i^^ 2 

i^' = ^ = -^{^ + V^ts'], (13) 

dr s 1 — 77, 



dt dr [s'y I — n 1 — n 1 — n 1 — n 

(14) 
With these substitutions we obtain from equation (1) the following reduced 
form of the Lane-Emden equation 

u\t) = ^u^ + {^-^^t + nu^ n^l. (15) 

n — 1 (n — 1)^ 

Where u{t) = s'{t) and u'{t) = s"{t). This highly nonlinear first order 
differential equation has solutions for n = and 5. These correspond to 
the well known solutions for 'ip namely — (r^/6) + (l/6x) and (1 + r^/3)~^/^ 
respectively. The only other known solution to (15) is the spherical Bessel 
function corresponding to n = 1 [5], [6] . To investigate the solution for other 



permissible values of n we first transform the variable u{t) as —l/y{t). This 
substitution weakens the nonlinearity of the equation as 

ht + 1" 
y'{t) = a-—^, n^l, (16) 

where a = (n — 5)/(n — 1) and h = 2(3 — n)/{n — 1)^. For the infinitesimals 
(^, 0)The symmetry condition for equation (16) gives 



6t + r, , bt + t"^,^ , b + nr~^, ,M + t 



0„+(0„-^„)(a-^^)-^„(a-^^) =^( ) + ,, ^2 ^ 

(17) 
For the infinitesimals of this equation we find that no non-trivial solution to 
(16) exists for the above ansatz. 

However one immediate consequence of the infinitesimals in (9) is an 
invariant particular solution for the Lane-Emden equation. We first have 



t, n — 1 r [n 



*'W = :^ = z^(^). *"M = 9^^^^^ (18) 



Substitution into equation (1) we obtain the particular solution in an implicit 
form 



r = ;^^^fj. "^1- (19) 



2{n-3)ilj 
(n — 1)2 r^ 

For n less than one this solution is bounded for all values of r. For n >~ 1 the 
solution has a singularity at r = 0, which becomes stronger as n increases. 
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